Rotational stacking and its electronic effects on graphene films grown on 4H-SiC(0001) 
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We examine the stacking order of multilayer graphene grown on the SiC (0001) surface using low- 
energy electron diffraction and surface X-ray diffraction. We show that the films contain a high 
density of rotational stacking faults caused by three types of rotated graphene: sheets rotated 30° 
and ±2.20° relative to the SiC substrate. These angles are unique because they correspond to 
commensurate phases of layered graphene, both with itself and with the SiC substrate. Ab intio 
calculations show that these rotational phases electronically decouple adjacent graphene layers. The 
band structure from graphene at fault boundaries displays linear energy dispersion at the A'-point 
(Dirac cones), nearly identical to that of a single graphene sheet. 
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In the last few years an intriguing series of experi- 
ments suggests that a new all-carbon paradigm for elec- 
tronic circuits may be possible 0, 0, 0,13] ■ In this system 
graphene sheets are lithographically cut into ribbons to 
produce gates and wires from a single material [H, [H, @ . 
The promise of this new approach to electronics rests on 
the ability to make large single or multilayer graphene 
sheets on a substrate that preserves the electronic prop- 
erties of an isolated graphene sheet. 

Both exfoliated graphene flakes 0, 0] and multilayer 
graphene grown on SiC 0, S] exhibit 2D transport prop- 
erties characteristic of chiral Dirac electrons expected 
for an isolated graphene sheet. These include a Berry 
phase of tt in the integer quantum Hall effect and sup- 
pressed back-scattering 0, 0, E|. SiC-grown films of- 
fer the most practical and scalable approach to 2D 
graphene electronics, but the effect of the substrate 
and graphene stacking must be clarified. Recent X-ray 
diffraction measurements and electronic structure cal- 
culations showed that the graphene-SiC(OOOl) substrate 
interaction is weak, except for a strongly-bound (non- 
conducting) buffer layer isolating subsequent graphene 



layers from the SiC [10L lllj . However, the remaining 



paradox is that normal AB. . . stacked (Bernal) graphite 
breaks the equivalency of A and B atoms in a graphene 
sheet [12], [l3| so that multilayer films should not exhibit 
the graphene-like properties that are clearly observed on 
4H-SiC(0001) substrates HQ. As we wil1 show : in this 
case nature provides a new stacking sequence that pre- 
serves the electronic symmetry of graphene. 

In this Letter we present Low Energy Electron Diffrac- 
tion (LEED) and surface X-ray scattering data for mul- 
tilayer graphene grown on the 4H-SiC(000f) carbon- 
terminated surface (C-face). We show that graphene 
grows in three forms on this surface: layers rotated 30° 



(#30), or ±2.20° {R2 ± ) from the SiC bulk [1010] direc- 
tion. In contrast, 4H-SiC(000I) (Si-face) films orient only 
in the i?30 phase (known as the (6\/3 x 6^3)^30° recon- 
struction in SiC coordinates) [7J]. X-ray diffraction on 
C-face films confirms that all three rotated phases occur, 
causing a high density of rotational fault boundaries be- 
tween i?30 and R2 ± layers. Ab initio electronic calcula- 
tions for this type of stacking show that adjacent rotated 
planes become electronically decoupled, preserving the 
Dirac dispersion at the if -point. These results explain 
how, even in multilayer graphene, the films maintain the 
electronic properties of an isolated graphene sheet. 

All substrates were 4H-SiC(0001) prepared as previ- 
ously reported [§]. Once samples are graphitized, they 
remain relatively inert, allowing them to be transported 
into the vacuum chamber. The X-ray scattering experi- 
ments were performed at the Advanced Photon Source, 
Argonne National Laboratory, on the 6IDB and C-/1CAT 
beam lines at 16.2 keV photon energy. Samples were 
mounted in a vacuum cryostat in the diffractometer. Re- 
ciprocal space points are reported in the reciprocal lat- 
tice units (r.l.u.) of the standard graphite hexagonal re- 
ciprocal lattice, q= (ha* ,kh* ,£c*), where |a*| = |b*| = 
27r/(av / 3/2) and |c*| — 2tt/c). The nominal lattice con- 
stants for graphite are a = 2.4589A, c = 6.674A [l5|. 

While it is known that graphene grows epitaxially only 
in the R30 phase on the 4H-SiC(000I) Si-face, multilayer 
graphene grown on the C-face was thought to have a 
high degree of azimuthal disorder because of streaking in 
LEED images [7]. However, a detailed look at the diffrac- 
tion shows that the rotational disorder is not random. 
This is demonstrated in the LEED image from a film 
with ~ 10 graphene layers [Fig. HJa)]. The pattern has 
two characteristics; (i) an oriented i?30 film evidenced by 
graphene spots (rods) rotated ±30° from the SiC [1010] 
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FIG. 1: (a) A LEED image acquired at 67.9eV from 4H- 
SiC(OOOl) C-face with ~ 10 graphene layers showing only 
graphite spots and diffuse rings. The SiC [1010] direction 
is marked for reference. The SiC (6v3 x 6y/S)RS0° unit vec- 
tors A* and B* are shown for reference, (b) X-ray azimuthal 
scans of the diffuse graphite ring around = 0.0 and \q\ = |a*|. 



direction and (ii) azimuthally diffuse rings centered at 0° 
from the SiC [1010] direction. Note that the diffuse rings 
are not continuous but split. This is seen more clearly in 
Fig. CJb) that shows an X-ray azimuthal scan (<j> scan) 
taken at the radial position of the graphite rod (\q\ = a*) 
in Fig. nja)) around = 0.0°. The scan shows diffuse 
intensity that is peaked at <j) = ±2.2°. 

The significance of the ±2.2° preferred rotation is two- 
fold. First, graphene is nearly commensurate with the 
SiC substrate on a length scale equal to a (13 x 13) 
graphene unit cell (rotated 30° from SiC). The (13x13) 
cell is ~ 0.14% smaller than the SiC (6^/3x6^)^30° 
cell. What has not been recognized is that there are 
two additional ways to orient a (13 x 13) graphene sheet 
that give rise to a film commensurate with the SiC 
(6v / 3x6v / 3)-R30° structure. They can be calculated when 
the mag nitude of the SiC (6^3 x 6%/3) reciprocal lattice 
vector is nearly equal to the graphene reciprocal lattice; 



\nA* + mB* 



(1) 



n and m are integers. A* and B* are the reciprocal lattice 
vectors of the (6voX 6\/3)-R30° structure (see Fig.QJa)), 



\A* 



\B* 



1/(6^3) with a* SiC = 2.3554A .The 



rotation angle of the commensurate graphene relative to 
SiC can be calculated for different m and n's: 
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FIG. 2: Schematic (V13 x V13)-R46.1 fault pair unit cell 
(dashed line). Dark circles are R30 graphene atoms (a and b 
are graphene unit vectors). Gray circles are graphene atoms 
in the R2 + plane below, rotated 32.204° from the top plane. 



Equation ((T|) is satisfied when (n, m) — (13, 0), (8,7) or 
(7,8). All three solutions give the SiC (6^3x6^3)^30° 
cell but the first solution has graphene rotated 30° rela- 
tive to SiC while the other two solutions have graphene 
rotated ±2.204° relative to SiC. Graphene grown on the 
Si-face of SiC only locks into the 30° structure. How- 
ever, as the diffraction in Fig.[lja) clearly shows, all three 
phases appear on C-face grown multilayer graphene. For 
future discussion we will index the two spots near </> = 
as the (1,0,£) R2 + an d (1,0,^)^2- graphene rods. It is 
worth noting that, aside from the (1,0,^)^ 2 ± rods, other 
(6-^/3x6-^/3)^30° reconstruction spots must be weak since 
they are not seen in LEED. This is in contrast to Si-face 
grown graphene where the additional spots indicate a 
complicated interfacial graphene reconstruction. Appar- 
ently C-face graphene has a weaker substrate interaction 
compared to Si-face graphene that may explain why the 
additional rotated phases are specific to C-face films. 

The significance of these three phases is even more im- 
portant if we recognize that two graphene sheets can 
be rotated relative to each other in a number of ways 
that make the two sheets commensurate with each other 
The lowest energy commensurate rotation angles 
are precisely <£ = cos _1 (ll/13) or cos" 1 (-11/13)- 120°) 
i.e, 30 ± 2.204° 16]. This bi-layer commensurate struc- 
ture corresponds to a graphene (vT3x ^/l3)R(±46.1°) 
cell that costs 3-5meV/atom more than AB. . . stacking 
[a schematic of a fault pair is shown in Fig. [2] [l6( • 

While the observation of three rotational phases is in- 
teresting, it is the stacking of these rotated planes that 
bears directly on their electronic properties. The i?30 
and R2 ± phases do not exist as isolated domains. Instead 
all three rotations are present in a multilayer graphene 
stack and lead to a high density of R30/R2 ± fault pairs. 
The most direct evidence for this comes from high- 
resolution X-ray diffraction. Figure [21a) shows X-ray ra- 
dial scans [see Fig. [Ha)] through the graphite (1,0, £) 
rod for different values of i. The two non-dispersing 
peaks correspond to a normal graphene (1,0, £) surface 



rod and the other to a surface rod (1 + Ah,0,£) from 
graphene with a compressed in-plane lattice constant. 
The peak separation corresponds to a lattice compres- 
sion of Aa/a=Ah/(l + Ah) = -0.28 ± 0.01%. Note that 
the compressed and uncompressed rods widths are the 
same, meaning that the ordered size of these two types 
of graphene sheets are similar. It is important to real- 
ize that the compressed lattice is not a result of epitaxial 
strain between the graphene and the SiC substrate, which 
would increase the graphene in-plane lattice constant. 

We can identify the compressed graphene as those 
sheets at the R30/R2 ± fault boundary shown in Fig. [2] 
This conclusion is derived from two key pieces of infor- 
mation found in the intensity modulation of the (1,0, £) 
rod [see Fig.[3Jb)]. First, in normal graphite, the primary 
fault type is rhombohedral ABC. . . stacking. Faults of 
this type would produce a peak in the (1, 0, £) at £ = 1.5 
that is clearly not seen in the data in Fig.^b). 

Second, A A. . . and ABC. . . faults are expected to 
cause small inter-layer contractions ~ 0.2% (l7j . Instead, 
our X-ray results show a large interplanar expansion at 
the fault boundary. This is shown in the inset of Fig. E2b) 
where the experimental peak at £ = 2 is shifted to a 
slightly lower value. This indicates a larger inter-layer 
spacing than bulk graphite. We can estimate the inter- 
planar expansion at the fault using a model where ran- 
dom rotational faults are introduced with a probability 
7, and each R30/R2 ± boundary expands by e. A fit with 
7 = 0.38 and e = 0.06A is shown in Fig.^b). This large 
expansion (1.8%) is characteristic of azimuthally disor- 
dered turbostratic graphite with many rotational faults 
that cause significant interference of n* states between 
rotated planes [la ]. Note that this large interlayer ex- 
pansion is also consistent with the in-plane contraction 
of the fault. Graphite's thermal expansion is negative in- 
plane and positive out-of-plane [l8j] . A weaker inter-layer 
bond caused by the large expansion at the fault allows 
the in-plane bonds to contract [3, [n| • 

We expect significant changes in the electronic prop- 
erties of the graphene films when the rotated stacking 
fault density is as high as we observe. For a R30/R2 ± 
fault pair there are only 2 atoms/sheet out of 52 in the 
(\/l3 x \/l3)i?(±46.1 ) cell that are in high symmetry 
positions. This suggests weak interplanar interactions in 
the fault pairs. To understand how the rotational stack- 
ing affects the electronic properties of these films, we have 
performed an ah initio Density Functional Theory (DFT) 
band structure calculation for a R30/R2 ± free standing 
graphene fault pair. The calculations are performed us- 
ing the VASP 1201] code within the generalized grad ient 
approximation [2lJ . Ultra soft pseudopotentials [22j are 
used with a plane wave basis cutoff equal to 211 eV. All 
calculations are performed on the same bi-layer commen- 
surate cell, the (v / 13x\/l3)i?(±46.1 ). This cell contains 
two graphene sheets of 26 carbon atoms each, rotated 
32.204° relative to each other [see Fig. [2]. The empty 



0.05 




FIG. 3: (a) Radial (h) scans through the graphite (h = 1, 0, £) 
rod for different I (see Fig. (TJ. Scans show two peaks cor- 
responding to the normal graphene (1,0, •£) rod and a com- 
pressed graphene (1 + Ah, 0, 1) rod. (b) Integrated intensity 
of the (1,0, 1) rod. Solid line is a fit to a random rotational 
fault model. Vertical bar marks the expected position of a 
peak from ABC. . . faults. Inset is an expanded view near 
£ = 2 showing a shift to smaller I (larger inter-layer spacing). 



space width is equal to 24 A. The total energy of the ro- 
tated stacking cell is 1.6meV/atom higher than an AB. . 
cell. This energy difference is slightly smaller than other 
estimates of this structure [lj| . 

For bilayers, the interlayer distance is fixed to 3.39A. 
As a check, we have varied this distance and found no 
qualitative change in the band structure results. This 
check was performed since DFT is known to poorly de- 
scribe van der Waals forces and gives rise to theoreti- 
cal graphene interlayer spacings significantly larger than 
experimental values. However, we point out that the 
C-short ultrasoft pseudopotential used here has been ex- 
tensively tested [ill HH and was shown to correctly de- 
scribe the band structure of graphite . Integration over 
the Brillouin zone is performed on a 30 x 30 x 1 grid 
in the Monckhorst-Pack scheme to ensure convergence 
of the Kohn-Sham eigenvalues. Due to the faulted cell 
symmetry, the if -points of the graphene lxl cell for 
the 2 layers of the rotated stacking are translated to the 
If -point of the (\/l3x \/l3)i?(±46.1 ) Brillouin zone. 

The results of this calculation are shown in Fig. [D 
where we compare the band structure for an isolated 
graphene sheet, a graphene bi-layer with Bernal AB. . . 
stacking, and a bilayer rotation fault pair with the 
(v / 13x v / T3)i?(±46.1°) structure (the T K M direction 
shown is the (\/T3xv / T3)i?±46.1 cell high symmetry di- 




FIG. 4: Calculated band structure for three forms of 
graphene. (i) isolated graphene sheet (solid line), (ii) AB. . . 
graphene bilayer (dashed line), and (iii) R30/R2 + fault pair 
(dots). Inset shows details of band structure at the if-point 
showing no difference between the Dirac cone for an R30/R2 + 
fault pair and a single graphene layer. 



rcction) . The main differences in the electronic structure 
among the three forms of graphene shows up in the dis- 
persion curves in the vicinity of the if-points. The band 
structure for an isolated graphene sheet shows the known 
linear gapless dispersion (Dirac cone) of the it bands at 
the JsT-point. The normal Bernal AB. . . stacking of 
graphene breaks the sublattice symmetry, giving rise to 
splitting of the 7r bands with a corresponding change to 
a parabolic shape and a lower group velocity [13| . With 
the rotated fault, the linear dispersion is recovered in 
the vicinity of the i^-points. This dispersion is identi- 
cal to the graphene dispersion (same Fermi velocity) and 
clearly shows that in the rotated layers, the atoms in 
the A and B sublattices are identical. This result also 
holds for infinite stacks: a graphite-like system made of 
graphene sheets rotated alternately by 0° and 30± 2.204° 
also shows a linear dispersion near the if-point. This re- 
sult is similar to the continuum description of Santos et. 
al 24j for two graphene sheets rotated by much smaller 
relative angles. 

In conclusion, we have shown that multilayer graphene 
grown on the carbon terminated face of 4H-SiC does 
not grow as a simple AB. . . stacked film. Instead the 
graphene-SiC and graphene-graphene commensurations 
produce a high density of rotational faults where ad- 
jacent sheets are rotated 30° ± 2.204° relative to each 
other. These rotational faults cause adjacent graphene 
sheets to decouple electronically. The result is that the 
band structure of faulted sheets is nearly identical to iso- 
lated graphene. Specifically, the Dirac dispersion at the 
isf-point is preserved even though the film is composed 
of many graphene sheets. This may explain why mag- 
netotransport [4j and infrared magnetotransmission 14 1 
experiments on similar samples give results very similar 
to those of an isolated graphene sheet. 
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